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Abstract 



We investigate the interface dynamic in Laplacian growth model, using the 
conformal mapping technique. Starting from the governig equation for the 
conformal map, obtained by B.Shraiman and D.Bensimon, we derive different 
possible forms of the equation. Some of them have evident physical represen- 
tation and are convenient for computer simulation, the results of which are 
presented in this paper. 

Formation and evolution of interfaces in systems with phase transition of the first kind 
led to one of the most striking phenomena of self-organization - the growth of dendritic 
structures. These phenomena are inherent to many branches of physics, chemistry and 
biology and are quite similar from mathematical point of view. They are known as the 
moving boundary problems or Stefan problems 

Such problems, as usual, have no analytical solutions caused by the essential non-linearity 
stipulated by the representation of the boundary conditions for parameters, allocated in 
space on the free boundary. In the given class of systems with free boundaries it is possible 
to single out a subclass where the distribution of the field in the space is stationary. It 
is known as quasistationary Stefan problem. Usage of the methods of conformal mapping 
allows to reduce the last problem to a dynamics of some physical fields in the system with 
fixed boundaries. This approach is more suitable for analytical and computer simulation. 

Theoretical methods of the functions with the complex variable appear to be rather 
useful in the construction of the effective surface dynamics of the interface at least in the 
two-dimensional problem. In this case the central part belongs to the conformal mappings 
of the physical regions of the space onto some special kind of regions. A large number of 
models for interfacial dynamics have been proposed by now. Reviews of these problems 
can be found in [^j-^. The most famous equation has been introduced by Shraiman and 
Bensimon and now is widely known as the conformal mapping equation and it is still 
under investigation. 
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Let us consider the problem of Laplacian growth, which can be formulated in the following 
way. We investigate the growth of two-dimensional region D and suppose, that its boundary 
dD, which represents the physical interface, is the analytic Jordan curve F 0. The field Lp 
outside of the region D, obeys the Laplace equation 

Aif = 0. (1) 

The boundary dD grows at a rate, that is proportional to the normal gradient of the field of 
the interface. Therefore the evolution of the interface is governed by the following equation 

Vn = -VnV^lr ■ (2) 
The potential field (p satisfies the following boundary condition at the interface 

(/?|p = do^lr > (3) 

where k is the curvature of the interface F, do - is the dimensionless surface tension param- 
eter. 
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FIGURES 

FIG. 1. Time-dependent conformal mapping from the exterior of the unit disc onto the domain 
of our interest. 



Equations (|T])-® determine the boundary- value problem with moving interface, or Stefan 
problem. Usage of the the methods of conformal mapping @-12,14|, which are based on the 
Riemann mapping theorem, allows to simplify the considered task. The theorem states the 
existence of a conformal map from the exterior of D onto a standard domain, for example, 
the exterior of the unit disk, and boundary dD corresponds to unit circle. In outline one 
can parametrize the evolution of the interface in z plane by the time-dependent conformal 
map F {w,t), which takes the exterior of the unit disk at each instant of time, \w\ = 1, onto 
the exterior of D. Therefore the evolution of the interface can be represented in the terms 
of 



F {e, t) = lim F{w,t), 0<e <27v 



(4) 



As it is shown by Shraiman and Bensimon [§, in case of free surface tension {do = 0) the 
conformal map fulfills the following equation of motion 

dT{e,t) 



dT{e,t) _ . dT{e,t) ^ 



dt 



de 



de 



(5) 



where S [■] stands for Schwartz operator. For the exterior of the unit circle with some 
boundary condition / {6) the Schwartz operator can be represented in the following way 



2tt 



S[fm = -^ld9 





, e + w 



/ {0') + zC, 



(6) 



dT(e,t) 



where C is an arbitrary constant. Denoting the action of the Schwartz operator on 
as the limit of an analytic function G{w,t), which explicit form is unique and explicitly 
determinable the equation of motion for F {w, t) (^) is transformed to 

dF{w,t) dF{w,t) 



w- 



-G{w,t). 



dt dw 

Let us also write down the logarithmic derivative of the with respect to w 



(7) 



d ^^dF{w,t) _ fdF{w,ty 
dt dw \ dw 



d ( dF{w,t) ' 

— w — G{w,t) 

dw \ dw 



(8) 



Using the explicit forms of the conformal map F {w,t) one can reduce the partial differential 
equation governing the interface dynamics, to a set of ordinary differential equations for 
the dynamics of the critical points of the conformal map This facilitates the analytical 
and numerical investigations of the complex behavior in the evolution of the interface Ij^-Q . 

In spite of this fact sometimes it is convenient to introduce the following field |]5|, p!TlJT2 
for these investigations 

dF{w,t) 



rj{e,t) 



dw 



(9) 



which has clear representation - the stretch of the interface F {6, t) under conformal map 
F^^ {z, t). The conformality of F {w, t) demands that the critical points of the map, zeroes 
of its derivative, all lie within the unit disk. Therefore In '^'^j^'*-' is an analytic function 
outside of the region D. Regarding that 



In 



dF {w,t) 
dw 



In 



dF {w,t) 



dw 



(10) 



and solving Schwartz problem outside the unit disc in w plane, one can recover the equation 
of motion for F {w,t) from the equation of motion for i] {0,t), as well as from ip {6,t) one 
('0(6',t) represents oriented angle between tangent vectors Tw and tf(^ui) P- To keep clear 
the sense of using (H) we should also write down the relations between rj {0,t) and tfj {9,t) 
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iP(e,t) = [ d9' coth 
2tt J 



e'-e 



\nr]{e',t) 



(11) 



In the present paper we investigate the properties of the equation of motion, obtained 
by Shraiman and Bensimon (|^), and derive the equation of motion for r] {0,t), that is more 
suitable for numerical simulation. In this case let us note that 

-2 



dr {e,t) 


~2 


dF {w, t) dw 


-2 


dF{w,t) 


de 




dw dO 




dw 



(12) 



thus G{w,t) = S {t]'^ {9 , t)) , where all constants in Schwartz operator are explicitly deter 
minable . Regarding the real part of the left side of the equation (H) as ^ In 
and simplifying the right side of the equation one can get 



dF(w,t) 
dw 



d_ 
It 



In 



dF{w,t) 



dw 



Re 



^. , dG(w,t) ^. d ^ dF(w,t) 

G (w, t) + w + wG (w, t) — In 

dw dw dw 



(13) 



Determining the real part for each term within square brackets of the right side of ([131) for 

if) 

w = e 



Re [G{w,t)]U^.e = v' iO,t) 



Re 



w- 



d G{w,t) 
dw 



= ^Im[G(«;, 



(14) 
(15) 



Re 
Re 



wG (w, t) —— In ; 

dw dw 



d ( 

wG {w, t) — In 
dw \ 



dF{w,t) 



dw 



-v'{0,t)^^l^ + lm[Giw,t)]-[^ln 



+ iip {w, t) 

dF {w,t) 



dw 



(16) 



and taking into account (|^), (0) we rewrite the equation (|T3|) in the following form 



d 



de 



+ Im [G{wM^=e^e-\nr'{e,t). 



(17) 



Taking into account (|TT]) and the fact, that imaginary part of Schwartz problem solution 
can be written in the terms of 



Im [G {w 



5 yj \'w=e^' 



^'j de' coih^-^rf {9\t) , 





(18) 



one can simplify equation (0) and obtain the following equation of motion for rj {6, t) 

= -r,^ (6, t) - rj (9, t) ^i- / d9' coth ^r;^ (9', t) + 





, Oil (g, t) 1 J 9' -9 2 

+^^^jd9 coth^r, i9,t). 





(19) 



The last equation has similar form to the derived one, using the variational principles of 
conformal maps P JTl| , p!2[] . 



In order to derive the proper equation of motion for i] {6, t) in the presence of surface 
tension {do 7^ 0) we should notice, that all we need is to obtain the proper form of G {w, t) 
for the nonzero surface tension [|T^]. 

Thus, one can rewrite the desired evolution equation for mapping (|^) in the following 
form: 



dF{w,t) dF iw,t) ^ 

—W ; Go [W, t) 



dt 



dw 



(20) 



where Gq {w, t) is an analytic function of w, which real part on = 1 is specified as follows 



Re [Go {w,t)] 



dw 



W 



dF{w,t) 
dw 



(21) 



Comparing equations (|20|) and (0) we should notice, that to derive the equation of mo- 
tion for the properly introduced 77 {9, t) , we have to define the real and imaginary parts of 
Go 



Re [Go [w, t)] L^,.e = {0, t) (l-do§^ Im [S [k {w, t)]] \ 



(22) 



27r 



Im [Go {w 



where 



5 ^Ji \w=e'^'^ 



9' -6 



d 



_ / d9' coth ^^7^' (9', t)[l- doj^ Im [S [k (w, t)]] L^,,^, ) , (23) 



d9' 



2tt 



Im [S [k {w, t)]] l^^.e = ^ Jd9' coth ^-^k {9', t) 





(24) 



After the similar simplifications for the free-tension case one could obtain the following 
equation of motion for 77 {9, t): 



a 1 '1 , / .91" 



\ 



27V \ 27V 

+rf t) 1 1 - rfo A i_ I dO' coth (^', tU^^ JdO' coth ^ In (6', t) + 



/ 



de 27: J 2 ' ^ ' M ^de'2Ti J 2 ^ ' ^ I ' v ; 

\ / 

where /c {6, t) is the curvature of V {0, t) : 

k (9, t) = -r^ (e, t) = ^ ^) / ^ 1^ ^ (^'' ^) • (26) 



One can also notice that 



Vn{e, t)=r^{e,t){l-do^ Im [S [k {w, t)]] L^^.. ) (27) 



Equations (p5D-(p7D represent the system of intrgrodifferential equations that fully de- 
scribes the evolution of the interface. 

To solve the derived singular integral equations with Hilbert kernel, the useful approach 
is developed, based on interpolation of the expressions of integration in singular integrals 
by trigonometric polynomial. For example, for each number = 2n of the nodes we can 
use the formulae [|13| 



Ue) = ^ E m) sm[n{9 - e,)) cot ^ (28) 
27r .^^ 2 



Oj = ^, j=0,l,...,2n-l (29) 



which lead to the next integral representation: 



±- J de' coth ^/(«,) = ^ E m) cot (30) 



" "'^ i=o 



where 

Q 27rm + 1 



«m = — n ' m = 0, 1, ...,2r2 - 1 (31) 

Zn 
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The system of nonlinear integrodifferential equations was investigated numerically. The 
basis of the numeric simulation was the approximation of unknown quantities of functions 
by above mentioned trigonometrical polynomial. Expansions of this type make it possible to 
use well-known quadrature formulae to evaluate singular Integrals . As a result, the system 
of equations and ( [27| ) reduces to a system of ordinary non-linear equations with 

full Jacoblan. 

It is possible to reproduce geometrically defined interfaces in parametric form by using 
the equations: 



In numerical experiment the evolution of the form of an interface of phases in an assoti- 
ation from an aspect of initial perturbations 6r]\t=o of a surface energy do was investigated. 

Figures (|^) and @ illustrate the characteristic form of non-uniform distributions of 
the r], i/j and the corresponding interfaces for different initial perturbations. The reader 
should observe that relatively smooth variations of the interface correspond to quite marked 
variations in the field r] and hence also in the angle ip- Thus, relatively small errors in 
determining do not cause marked changes in the interface geometry. In addition, as is evident 
from the figures, relatively even long sections of the Interface contract when mapped into 
the interval [0,27r], whereas pronounced changes in the contour F, conversely, are stretched 
out. In this sense such mappings of onto are "adaptive" with respect to information about 
the detailed geometic structure of the interface. 

FIG. 2. Non-uniform distributions of the r/ {9,t) , tp id,t) and the corresponding interfaces for 
the following initial perturbation: r]{e,tQ = 0) = 0.1 - 0.01 cos(46l); do = 0.1; h = 500; ^2 = 1000; 
ts = 3000; t4 = 6000. Bold interfaces on the figure correspond to mentioned above instances of 
time. 




(32) 
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FIG. 3. Non-uniform distributions of the r] {0,t) , tp {d,t) and the corresponding interfaces for 
the following initial perturbation: ri{e,to = 0) = 0.1 - 0.01 cos(66'); do = 0.1; h = 500; t2 = 1000; 
^3 = 3000; ^4 = 6000. Bold interfaces on the figure correspond to mentioned above instances of 
time. 
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